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We derive the Onsager response matrix of fluctuation-mediated spin-collinear transport through a ferromag¬ 
netic insulator and normal metal interface driven by a temperature difference, spin accumulation, or magnetic 
field. We predict magnon-squeezing spin currents, magnetic field-induced cooling (magnon Peltier effect), tem¬ 
perature induced magnetization (thermal magnetic field) as well as universal spin Seebeck/Peltier coefficients. 


Finite temperature effects on the transport properties of 
magnetic nanostructures in attracts considerable attention 
since the discovery of the spin Seebeck effect (2^4) that 
thwarts conventional thermoelectrics. Of special interest are 
heterostructures of magnetic insulators such as yttrium iron 
garnets (YIG) with heavy normal metals such as Pt, where the 
latter, via the inverse spin Hall effect, function as spin cur¬ 
rent detectors. Here we report a linear response approach to 
thermal transport through interfaces between ferro- or ferri- 
magnetic insulators (FI) and normal metals (N) that extends 
our treatment of the spin Seebeck effect a to the spin Peltier 
effect and leads to the prediction of, e.g., a magnon Peltier 
effect and its Onsager reciprocal, a thermal effective field. 

The Fandau-Fifshitz-Gilbert equation for the dynamics of 
a magnetization in an effective magnetic field B 0 

m = — 7 m x (Bo + b) + am x m ( 1 ) 

is based on the assumption that the modulus of the spa- 
tiotemporal magnetization texture M(r, t) is constant, i.e. 
M(r,f) = M s m(r, t) and |m| = 1, which is valid at tem¬ 
peratures sufficiently below that of the magnetic phase tran¬ 
sition. The FFG predicts a temperature-induced reduction of 
the time-averaged equilibrium magnetization by considering a 
stochastic magnetic field b(r, t) that induces thermal fluctua¬ 
tions of m around the equilibrium direction. Thermal noise 
is characterized by the spatiotemporal correlation function 
( 6 j(r, t)bj(r', t')) that by the Fluctuation-Dissipation Theo¬ 
rem (FDT) can be expressed in terms of an integral over the 
Bose-Einstein distribution function of the magnon excitations 
and proportional to the Gilbert damping constant a @. Mi¬ 
croscopically, the magnetization noise in insulating ferromag- 
nets is caused by the magnetoelastic interaction that couples 
and equilibrates the magnetic and elastic order parameters. At 
interfaces to metals, spin pumping induces an additional en¬ 
ergy and angular momentum dissipation that increases the ef¬ 
fective damping and the magnetic fluctuations 12H9). 

Spin accumulations in the normal metal at interfaces to 
ferromagnets with transverse spin polarization generate spin- 
transfer torques mm, while the longitudinal ones have at 
zero temperature no effect on the magnetization. One might 
therefore, naively, expect that the magnetization of the insula¬ 
tor in FI|N bilayers is inert without outside spin injection or 


non-collinear magnetic fields. However, spin collinear trans¬ 
port phenomena in FI|N systems exist at finite temperatures 
by the magnetic thermal fluctuations that allow a longitudi¬ 
nal spin accumulation in N to act on instantaneous transverse 
magnetization components. The ensemble/time average of the 
thus induced spin currents is polarized along the equilibrium 
magnetization direction. 

Perturbation on a system at thermal equilibrium generates a 
response or “current” that is proportional to the “force” when 
the latter is sufficiently weak (Ohm’s Faw). In the presence 
of multiple forces and currents cross-correlations exist, ther¬ 
moelectrics being a prominent example. The linear response 
of such a system is then described by a “conductance” matrix 
that relates forces and currents, which possesses a fundamen¬ 
tal symmetry referred to as Onsager reciprocity Ifl2l that is 
very useful in spintronics na. Here we establish the Onsager 
matrix for transport through a normal metal and a ferromag¬ 
netic insulator contact (N|FI) that is actuated by a spin accu¬ 
mulation in N, temperature difference over the interface, and 
(pulsed) external magnetic field. Each matrix element rep¬ 
resents a different physical phenomenon, of which the spin 
Seebeck effect is just one 0H3IS1- The extended Onsager 
matrix discussed in the following has already been implic¬ 
itly used (without details and with reference to the present 
work) in the analysis of the spin Peltier effect El and in 
modelling spin Seebeck generators fT7l . The N|FI system 
has recently been analyzed by Bender et al. na using a 
“Golden Rule” treatment of the interface exchange interaction 
including angle-dependent spin transfer torques and quantum 
effects. However, this study does not take into account the 
magnetic field component parallel to the magnetization that is 
central to the present work. We focus on symmetry conserv¬ 
ing perturbations, thereby disregarding deterministic trans¬ 
verse spin accumulations and spin-transfer torques, which is 
allowed as long as the systems is far from the threshold of 
charge current-induced magnetic self-oscillations or magneti¬ 
zation reversal. All elements of the response matrix are then 
scalars. Recently, Nakata et al. nsi derived the Onsager ma¬ 
trix for a bilayer of magnetic insulators actuated by inhomo¬ 
geneous magnetic fields and temperature differences. 

Fig -0 sketchs the ferromagnetic insulator (FI) with nor¬ 
mal metal (N) contact. The equilibrium FI magnetization 
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FIG. 1: (Color online) Spin and energy current driven by the thermal 
bias (AT = Tp — Tv), spin chemical potential V s , and an external 
magnetic field B app at an FI|N interface. 
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is parallel to the uniaxial anisotropy field B 0 = T 0 z = 
(<*> 0 / 7 ) 2 1| (m). We adopt the “three reservoir model” by as¬ 
suming that thermalization of spin waves in FI and electrons 
in N is sufficiently efficient. The steady state in the pres¬ 
ence of a temperature gradient is then characterized by the 
temperatures of phonon (FI and N), magnon (FI), and elec¬ 
tron (N) systems, thereby disregarding spin dependent elec¬ 
tron temperatures in N 1221 . Transport is generated by the 
differences in the thermodynamic variables on both sides of 
the interface. This situation is amenable to magnetoelectronic 
scattering theory of transport parameterized by the interface 
scattering matrix. This picture has a wider applicability in cir¬ 
cuit theories, in which the interfaces generate boundary condi¬ 
tions between “bulk regions” that may be described by quasi¬ 
equilibrium distribution functions l23l . For comparison with 
experiments, the present analysis is a (crucial) building block 
in simulating entire devices mm. 

For simplicity, we consider here the limit of small phonon 
(Kapitza) heat conductances, which allows us to discard the 
phonons altogether G3. The thermodynamic state of the in¬ 
sulating ferromagnet is then characterized by the temperature 
Tp only. The normal metal is at electron and phonon temper¬ 
ature T n . We include the option of having long-lived spin ac¬ 
cumulation, i.e., a chemical potential difference between the 
spin up and spin down electrons in the frame of the FI equi¬ 
librium magnetization with quantization axis along z. The 
thermodynamic forces are then the external magnetic field 
B a pp = -B ap pZ = (w app / 7 )z, longitudinal spin accumulation 
V s = V s z = (hjjj s / 2e)z (in units of volt) in N, and temper¬ 
ature difference AT = Tp — Tn across the interface. With 
T = (Tp + Tjv)/ 2, the rate of change of the free energy reads 

F = SV s + Of - ATT ap p, (2) 

where S (t) = J N s(r, t)dV = Si is the total spin (in units 
of electric charge |e|) in N with spin density s(r ,t), Q is 
the heat current entering N, M z (t) = M s f FJ m z (r,t)dV is 
the total magnetization in FI with local magnetization com¬ 
ponent M s m z (r,t). The volume of the FI is V - Ad with 
interface area A and thickness d. F includes contributions 
from the electrons (first two terms) and the magnetization (last 
term). Since Onsager symmetry holds when the entropy gen¬ 
eration rate S = F/T equals the sum of the product of cur¬ 


rents and forces, we identify in Eq. @ V a , -A T/T, -Tapp 
as thermodynamic driving forces with spin current density 
j s = [H/2e)S/A, energy current density j'q = Q/A, and 
magnetization dynamics M z /jA, as the conjugated flows, re¬ 
spectively. The linear response coefficient matrix in 
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must then satisfy the Onsager reciprocal relations m, which 
demand that our response matrix is symmetric. The matrix 
element L s q parametrizes the spin Seebeck effect I51IT41 151. 
Other elements represent various physical effects; L ss is the 
longitudinal spin conductance, Lq s governs the spin Peltier 
effect, Lqq is the heat conductance, while the signficance of 
the new matrix elements LQ m and L m Q is discussed below. 

Macrospin model - We first treat the macrospin limit in 
which m = M/M with M = M s V constant in space and 
described by the Landau-Lifshitz-Gilbert equation including 
interface torques: 


rii = — 7 m x B e ff + am x m — c/w s m x (m x z). (4) 


where B e ff = Bq + B app + b includes the internal field 
Bo||z and the applied field B app ||z. The last term is the in¬ 
stantaneous spin-transfer torque acting on the magnetization 
by the longitudinal spin accumulation in N QS1Q31- The 
Gilbert damping parameter a is the sum of the intrinsic ao 
and a' = 7 hg r / ( 47 t M s d), the interface damping l25l . where 
g r is the real part of the spin mixing conductance per unit area 
of the FI|N interface lf23l . Thermal effects are modeled by the 
stochastic fields b = b 0 + b', which by the FDT, is associ¬ 
ated with the damping a = op + a ' 01261. Sufficiently far 
from the Curie temperature, only the transverse x and y com¬ 
ponents of m, b, and b' matter. To lowest order, the Fourier 
components of rn x . y satisfy m“ = J2j Xijl^j (b j = x,y), 
introducing the magnetic susceptibility 


w = -(1 + q 2 ) 1 (u> 0 - iotu -iu \ 

(u — Uq){uj — cOq) \ *<*> ujo — iauij 

with u// = ±ujq/(1 ± ia). According to the FDT ll27ll . the 
magnetic fluctuations for kpT hu >0 satisfy 
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where (• • •) denotes the thermal ensemble average. The cor¬ 
responding random fields are white: 
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The auto-correlation for the interface-induced random field b' 
has the form Eq. <[T[) but with a —► a' and Tp —>• Tv- 

We are interested in the DC (ensemble or time averaged) 
spin and energy currents across the FI|N interface. The trans¬ 
verse components average to zero, while the z-componcnt in 
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Eq. (i) has a bulk contribution driven by B app and interfacial 
contribution driven by (-B app , V s , AT): 


M z M . 

^=^aoWapp(mx(mxz )) z +3l, (8) 

where (• • -) z is the thermal average of the ^-component, 

M . , 

3 = — (a m x m — 7m x b 

gA 

-a' (w s + w app ) m x (m x z)) z ( 9 ) 

with rh « — 71 TLX B 0 . The associated energy current equation 
M 


jn = —- (a'rh ■ m — 7m • b' 

7 A 

-a' (w s + w app ) z • (m x rii)) 


( 10 ) 


follows from the interface contribution (terms proportional to 
a! and b') of the energy change rate: dE/dt = (d/dt)(— B 0 • 
M). At zero temperature, m = z and b' = rh = 0, 
and both and jg vanish, as expected. At finite temper¬ 
atures, and in spite of (m)||z, and Jq are finite because 
m and b' are correlated. The relevant equal-time correlators 
(TOjTOj), ( rriib'j ), ( rriiiTij ), and (to,; 6 ') can be derived from: 


(mi(t)mj{ 0)) 
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Plugging these into Eqs. (| 9 j [T 0 | , we arrive at the linear re¬ 
sponse relation 


( 



2a' k B T 
A wq 


1 

1 


vw 0 


w 0 




( 12 ) 


Magnon model - The macrospin model above holds only in 
the presence of strong applied magnetic fields or nanomagnets 
smaller than the exchange length. Otherwise the thermal spin 
wave excitations and magnetization texture m(r, t) (ug JJ in 
Fourier space) may not be disregarded. The internal magnetic 
field should then be augmented by the exchange interaction 
B 0 —> B 0 + (D/jH) V 2 m where D is the spin wave stiff¬ 
ness. The stochastic fields bo(r, t) and b'(r,i) then depend 
on position and b'(r,t) = b'(y, z, t)S(x) acts at the inter¬ 
face. After linearizing and Fourier transforming Eq. 0 into 
frequency and momentum space, m*’ u = where 

the magnetic susceptibility x k ’“ takes the form Eq. ([5) with 
wo —> wo + ( D/H)k 2 . According to the FDT, the equilibrium 
magnetization fluctuations satisfy HU: 
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where <5^, = < 5 (w — w')< 5 (k — k'). The correlations for b 0 
and b' (or b') can be inferred from Eq. ( fTJ) . The Planck dis¬ 
tribution regulates the frequency integral over the continuum 


of magnon density of states. We note that the magnetic field 
dependence of the spin Seebeck effect at room temperature 
indicates a magnon cut-off lower than k B Tp/h. I28d30l . 

The spin and energy currents across the FI|N interface are 
still given by Eqs. when substituting a! —» a' = a!d , 

b' — > b', and M —>■ M/d = M S A. Using Eq. & all equal- 
time-positon correlators in Eqs. «IM can be infered from: 


(TOi(0,0)mj(0,0)) = n F J 0 (x 0 )6 ij , (14a) 

(to,(0, 0 ) 76 ' (0)) = -a , ^(m l (0,0)TO J (0,0)), (14b) 


where n F = 7 ?i/M s A 3 is the total number of spins in the 
volume A 3 = ( 47 r D/k B T) 3 / 2 and A is the de Broglie ther¬ 
mal wave length for magnons. Ji(x 0 ) = 2^[xJtt(xo + 
x) l /(e Xo+x — l)da; with Xq = hcoo/k B T F . The expressions 
hold to leading order in a. In the classical limit Xq <C 1, 
J\ —> Z/ +:i / 2 with Z n the Zeta function. Using these correla¬ 
tors in Eqs. (|9|[T0|) leads to the central result of this paper: 
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(15) 

with /3 _1 = k B T F . The response matrix is symmetric as 
required by Onsager reciprocity and invertible, i.e. the forces 
are linear-independent as long as oq U 0 . 

Discussion - The spin Seebeck effect represented by L s q: 


j* a =3Z*^k B AT~ 


AT 0.14/zJ 
0.1 K nr 2 


(16) 


specifies the longitudinal spin current induced by the temper¬ 
ature difference AT EHEIEI- The numerical estimates here 
and in the following are for the Pt|YIG system at T = 300 
K with parameters given in Table [I] The inverse of the spin 
Seebeck effect is the spin Peltier effect given by Lq s : 


Jq = 3Zs — k B T —14 — 


5 A 3 


U s 1.3 x 10 5 J 
0.1 /rV m 2 • s 


(17) 


U drives an energy current that cools/heats the magnons lfl6l . 

The spin current driven by spin accumulations and external 
magnetic field 


a' ( 2e \ 

j s = 7/sm^app 4“ L SS U S = 2 Z3 H ( 7-^app 4” j. 1) ) 

(18) 

vanishes with temperature since we disregarded quantum fluc¬ 
tuations and the magnon chemical potential in the ferromag- 
net fT8l . L sm in Eq. ( fl8) quantifies the spin current induced 
by shifting the spin wave gap or “squeezing” the magnon dis¬ 
tribution function. The spin conductance L ss in Eq. (18) de¬ 
scribes the spin current injection by a collinear spin accumu¬ 
lation V s as observed recently by Cornelissen et al. f33j . The 
magnon accumulation rate induced by this spin injection is 
described by L ms . 
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Parameter 

Value 

Unit 

Reference 

M s 

1.4 x 10 5 

A/m 

ii 

D 

8.5 x 10~ 40 

Jm 2 

[35l|36] 

ao 

1(T 4 



9r 

10 19 

1/m 2 

03 

A 

1.6 

nm 


a' 

0.1 

nm 



TABLE I: Parameters for YIG and YIG|Pt interface. 


The energy current driven by the external held via Lq tii 


3 Q = 3Zi—k f 


5 A 3 


■T^B app ~ 


B app 7.4 x 10 7 J 
1 T m 2 • s 


(19) 


with z' = (15^7/2 ^ 9Z 2 / 2/Z :V2 )/2. is of the same order of 
magnitude as the electron-phononic thermal conductance for 
A1|A1 2 0 3 interface of ~ 2 x 10 8 W/(m 2 • K) l38l . 

We can estimate the time scale of the transient response to 
a suddenly switched-on magnetic held. According to the first 
matrix element in Eq. © with w app = 7 .Bapp 


(rh z ) 


2 jh y 

M a A 3 ^ 0 Bpp 


‘ll'lF Z 3 &QUJ app . 


(23) 


On the other hand, the value of m z at thermal equilibrium 
depends on the spin wave gap hui q: 


(m z (w 0 )) = 1 - (ml) - (to 2 ) ~1- n F J 0 (x 0 ), (24) 

that is shifted by w app , therefore 


reflects what we call a magnon Peltier effect. It can be ob¬ 
served by a temperature change of the ferromagnet generated 
by the applied magnetic held, analogous to the spin Peltier ef¬ 
fect caused by a spin accumulation Ha. The reciprocal to the 
magnon Peltier effect is the magnetization dynamics induced 
by a temperature gradient via L m q 


= 3Z S ke&T , (20) 

similar to applying a magnetic held. The magnitude of this 
thermal effective magnetic field can be estimated by setting 
M z = 0: 


Bat 


a' 3Z 5/2 k B AT 
a 2Z 3/2 "/h 


a' AT 
a' + «o d 0.1 K 


60 mT, (21) 


where the prefactor is of the order of unity as long as d <C 
a'/ao (~ 1 for YIG). 

The conventional Seebeck coefficient refers to the open cir¬ 
cuit thermopower voltage induced by a temperature differ¬ 
ence. We may analogously dehne a spin Seebeck coefficient 
Ss in terms of the spin thermopower, i.e. the spin accumula¬ 
tion generated by a temperature bias AT in a spin-open cir¬ 
cuit: Ss = (V s /AT)j Z=0 . This coefficient is not very rele¬ 
vant for metals that act as spin sinks such as Pt, but it charac¬ 
terizes the efficiency for weakly spin-dissipating metals such 
as Cu 03. The thermally induced spin accumulation at zero 
spin current reads: 


3^5/2 k B AT 
4^3/2 


Ss 



( 22 ) 


with z = 3Z 5/2 /4Z 3/ 2 ~ 0.385. For single (or multiple) 
parabolic magnon bands Ss — 33 /iV/K is a universal con¬ 
stant that does not depend on the spin wave stiffness, uniaxial 
anisotropy, and (in the considered regime) temperature. 

The spin-magnonic interfacial thermal conductance K m = 
{jQ/^T) j;= o = ( Lqq - L s qLq s L~})/T: 


Km 


, a! k 2 B T F 
z A 3 h 


~ 0.76 x 10 s 


W 

m 2 • K 


d(m z (u} 0 )) 

dm z — ^app 

acjr) 


t • / \ ^-^app 

= »fLu(x„) — , 


(25) 


where Li is the PolyLog function. The magnon relaxation 
time to the new equilibrium is therefore 


Sm z _ h Lii(a;o) ^ h 

(m z ) ak F T 2 Z 3 a\fkffjThid 0 2^3 ’ 


where the approximate expression is valid for Xq = 
huj 0 /k B T <C 1. For YIG with / = w 0 /2tt = 10 GHz at 
room temperature, r ~ 2 ns. Judging from the experimen¬ 
tal magnetic field dependence of the spin Hall effect at room 
temperature, k B T should be replaced by a magnon spectrum 
cut-off of 30 K lf28]430ll . leading to the estimate r ~ 7 ns. At 
even lower temperatures, k B T in Eq. ( |26l > should be replaced 
by hui 0 , leading to an upper estimate of 50 ns. 

The non-equilibrium magnetization M s 5m z in Eq. ( [25] ) can 
be interpreted as a non-equilibrium magnon accumulation and 
the magnetic field as its driving force. This interpretation pro¬ 
vides the link to theories of the electrically or thermally in¬ 
jected magnon Bose condensate in which the self-organized 
magnon chemical potential plays a crucial role mHHI. The 
weak magnon-phonon interaction reported recently by Cor- 
nelissen et al. . [331 is encouraging that a long lived magnon 
chemical potential and condensate can be electrically gener¬ 
ated. The pulsed magnetic field experiments suggested here 
do not require such a chemical potential and should yield use¬ 
ful insights into the magnon distribution function. 

In conclusion, we studied the magnetization dynamics cou¬ 
pled to spin and energy currents through FI|N interfaces at 
finite temperature as driven by collinear magnetic fields, spin 
accumulations, and/or a temperature bias. The response in this 
configuration vanishes with thermal fluctuations of the mag¬ 
netization. We establish the Onsager reciprocal relations for 
these response functions. The elements in the Onsager matrix 
are identified as the spin Seebeck effect, the spin Peltier effect, 
and previously overlooked ones such as the magnon Peltier ef¬ 
fect, effective thermal field, and magnon squeezing that have 
still to be observed experimentally. We identified a (nearly) 
universal spin Seebeck thermopower of 33 /iV/K. 
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